Transcriptome Analysis of Compensatory Growth and Meat Quality Alteration after Varied Restricted Feeding Conditions in Beef Cattle.

Compensatory growth (CG) is a physiological response that accelerates growth following a period of nutrient limitation, with the potential to improve growth efficiency and meat quality in cattle. However, the underlying molecular mechanisms remain poorly understood. In this study, 60 Huaxi cattle were divided into one ad libitum feeding (ALF) group and two restricted feeding groups (75% restricted, RF75; 50% restricted, RF50) undergoing a short-term restriction period followed by evaluation of CG. Detailed comparisons of growth performance during the experimental period, as well as carcass and meat quality traits, were conducted, complemented by a comprehensive transcriptome analysis of the longissimus dorsi muscle using differential expression analysis, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and weighted correlation network analysis (WGCNA). The results showed that irrespective of the restriction degree, the restricted animals exhibited CG, achieving final body weights comparable to the ALF group. Compensating animals showed differences in meat quality traits, such as pH, cooking loss, and fat content, compared to the ALF group. Transcriptomic analysis revealed 57 genes and 31 pathways differentially regulated during CG, covering immune response, acid-lipid metabolism, and protein synthesis. Notably, complement-coagulation-fibrinolytic system synergy was identified as potentially responsible for meat quality optimization in RF75. This study provides novel and valuable genetic insights into the regulatory mechanisms of CG in beef cattle.

Improving Genomic Predictions in Multi-Breed Cattle Populations: A Comparative Analysis of BayesR and GBLUP Models.

Numerous studies have shown that combining populations from similar or closely related genetic breeds improves the accuracy of genomic predictions (GP). Extensive experimentation with diverse Bayesian and genomic best linear unbiased prediction (GBLUP) models have been developed to explore multi-breed genomic selection (GS) in livestock, ultimately establishing them as successful approaches for predicting genomic estimated breeding value (GEBV). This study aimed to assess the effectiveness of using BayesR and GBLUP models with linkage disequilibrium (LD)-weighted genomic relationship matrices (GRMs) for genomic prediction in three different beef cattle breeds to identify the best approach for enhancing the accuracy of multi-breed genomic selection in beef cattle. Additionally, a comparison was conducted to evaluate the predictive precision of different marker densities and genetic correlations among the three breeds of beef cattle. The GRM between Yunling cattle (YL) and other breeds demonstrated modest affinity and highlighted a notable genetic concordance of 0.87 between Chinese Wagyu (WG) and Huaxi (HX) cattle. In the within-breed GS, BayesR demonstrated an advantage over GBLUP. The prediction accuracies for HX cattle using the BayesR model were 0.52 with BovineHD BeadChip data (HD) and 0.46 with whole-genome sequencing data (WGS). In comparison to the GBLUP model, the accuracy increased by 26.8% for HD data and 9.5% for WGS data. For WG and YL, BayesR doubled the within-breed prediction accuracy to 14.3% from 7.1%, outperforming GBLUP across both HD and WGS datasets. Moreover, analyzing multiple breeds using genomic selection showed that BayesR consistently outperformed GBLUP in terms of predictive accuracy, especially when using WGS. For instance, in a mixed reference population of HX and WG, BayesR achieved a significant accuracy of 0.53 using WGS for HX, which was a substantial enhancement over the accuracies obtained with GBLUP models. The research further highlights the benefit of including various breeds in the reference group, leading to enhanced accuracy in predictions and emphasizing the importance of comprehensive genomic selection methods. Our research findings indicate that BayesR exhibits superior performance compared to GBLUP in multi-breed genomic prediction accuracy, achieving a maximum improvement of 33.3%, especially in genetically diverse breeds. The improvement can be attributed to the effective utilization of higher single nucleotide polymorphism (SNP) marker density by BayesR, resulting in enhanced prediction accuracy. This evidence conclusively demonstrates the significant impact of BayesR on enhancing genomic predictions in diverse cattle populations, underscoring the crucial role of genetic relatedness in selection methodologies. In parallel, subsequent studies should focus on refining GRM and exploring alternative models for GP.

Genome-Wide Association Analysis Identifies the PMEL Gene Affecting Coat Color and Birth Weight in Simmental × Holstein.

Coat color and birth weight, as easily selected traits in cattle, play important roles in cattle breeding. Therefore, we carried out a genome-wide association study on birth weight and coat color to identify loci or potential linkage regions in 233 Simmental × Holstein crossbred beef cattle. The results revealed that nine SNPs were significantly associated with coat color (rs137169378, rs110022687, rs136002689, Hypotrichosis_PMel17, PMEL_1, rs134930689, rs383170073, rs109924971, and rs109146332), and these were in RNF41, ZC3H10, ERBB3, PMEL, and OR10A7 on BTA5. Interestingly, rs137169378, rs110022687, rs136002689, Hypotrichosis_PMel17, and PMEL_1 showed strong linkage disequilibrium (r2 > 0.8) and were significantly associated with coat color. Notably, Hypotrichosis_PMel17 and PMEL_1 were located in the gene PMEL (p = 2.22 × 10-18). Among the five significant SNPs associated with coat color, the birth weight of heterozygous individuals (AB) was greater than that of homozygous individuals (AA). Notably, the birth weight of heterozygous individuals with Hypotrichosis_PMel17 and PMEL_1 genotypes was significantly greater than that of homozygous individuals (0.01 < p < 0.05). Interestingly, the two loci were homozygous in black/white individuals and heterozygous in gray/white individuals, and the birth weight of heterozygous brown/white individuals (43.82 ± 5.25 kg) was greater than that of homozygous individuals (42.58 ± 3.09 kg). The birth weight of calves with the parental color (41.95 ± 3.53 kg) was significantly lower than that of calves with a non-parental color (43.54 ± 4.78 kg) (p < 0.05), and the birth weight of gray/white individuals (49.40 ± 7.11 kg) was the highest. Overall, PMEL appears to be a candidate gene affecting coat color in cattle, and coat color may have a selective effect on birth weight. This study provides a foundation for the breeding of beef cattle through GWAS for coat color and birth weight.

Genetic Origin and Introgression Pattern of Pingliang Red Cattle Revealed Using Genome-Wide SNP Analyses.

The Pingliang red cattle, an outstanding indigenous resource in China, possesses an exceptional breeding value attributed to its tender meat and superior marbling quality. Currently, research efforts have predominantly concentrated on exploring its maternal origin and conducting conventional phenotypic studies. However, there remains a lack of comprehensive understanding regarding its genetic basis. To address this gap, we conducted a thorough whole-genome analysis to investigate the population structure, phylogenetic relationships, and gene flows of this breed using genomic SNP chip data from 17 bovine breeds. The results demonstrate that Pingliang red cattle have evolved distinct genetic characteristics unique to this breed, clearly distinguishing it from other breeds. Based on the analysis of the population structure and phylogenetic tree, it can be classified as a hybrid lineage between Bos taurus and Bos indicus. Furthermore, Pingliang red cattle display a more prominent B. taurus pedigree in comparison with Jinnan, Qinchuan, Zaosheng, Nanyang, and Luxi cattle. Moreover, this study also revealed closer genetic proximity within the Chinese indigenous cattle breed, particularly Qinchuan cattle, which shares the longest identical by descent (IBD) fragment with Pingliang red cattle. Gene introgression analysis shows that Pingliang red cattle have undergone gene exchange with South Devon and Red Angus cattle from Europe. Admixture analysis revealed that the proportions of East Asian taurine and Chinese indicine in the ancestry of Pingliang red cattle are approximately 52.44% and 21.00%, respectively, while Eurasian taurine, European taurine, and Indian indicine account for approximately 17.55%, 7.27%, and 1.74%. Our findings unveil distinct genetic characteristics in Pingliang red cattle and attribute their origin to B. taurus and B. indicus ancestry, as well as contributions from Qinchuan cattle, South Devon, and Red Angus.

Genetic Association Analysis of Copy Number Variations for Meat Quality in Beef Cattle.

Meat quality is an economically important trait for global food production. Copy number variations (CNVs) have been previously implicated in elucidating the genetic basis of complex traits. In this article, we detected a total of 112,198 CNVs and 10,102 CNV regions (CNVRs) based on the Bovine HD SNP array. Next, we performed a CNV-based genome-wide association analysis (GWAS) of six meat quality traits and identified 12 significant CNV segments corresponding to eight candidate genes, including PCDH15, CSMD3, etc. Using region-based association analysis, we further identified six CNV segments relevant to meat quality in beef cattle. Among these, TRIM77 and TRIM64 within CNVR4 on BTA29 were detected as candidate genes for backfat thickness (BFT). Notably, we identified a 34 kb duplication for meat color (MC) which was supported by read-depth signals, and this duplication was embedded within the keratin gene family including KRT4, KRT78, and KRT79. Our findings will help to dissect the genetic architecture of meat quality traits from the aspects of CNVs, and subsequently improve the selection process in breeding programs.

Characterization of Chromatin Accessibility in Fetal Bovine Chondrocytes.

Despite significant advances of the bovine epigenome investigation, new evidence for the epigenetic basis of fetal cartilage development remains lacking. In this study, the chondrocytes were isolated from long bone tissues of bovine fetuses at 90 days. The Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) and transcriptome sequencing (RNA-seq) were used to characterize gene expression and chromatin accessibility profile in bovine chondrocytes. A total of 9686 open chromatin regions in bovine fetal chondrocytes were identified and 45% of the peaks were enriched in the promoter regions. Then, all peaks were annotated to the nearest gene for Gene Ontology (GO) and Kyoto Encylopaedia of Genes and Genomes (KEGG) analysis. Growth and development-related processes such as amide biosynthesis process (GO: 0043604) and translation regulation (GO: 006417) were enriched in the GO analysis. The KEGG analysis enriched endoplasmic reticulum protein processing signal pathway, TGF-ß signaling pathway and cell cycle pathway, which are closely related to protein synthesis and processing during cell proliferation. Active transcription factors (TFs) were enriched by ATAC-seq, and were fully verified with gene expression levels obtained by RNA-seq. Among the top50 TFs from footprint analysis, known or potential cartilage development-related transcription factors FOS, FOSL2 and NFY were found. Overall, our data provide a theoretical basis for further determining the regulatory mechanism of cartilage development in bovine.

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle.

BACKGROUND: A detailed understanding of genetic variants that affect beef merit helps maximize the efficiency of breeding for improved production merit in beef cattle. To prioritize the putative variants and genes, we ran a comprehensive genome-wide association studies (GWAS) analysis for 21 agronomic traits using imputed whole-genome variants in Simmental beef cattle. Then, we applied expression quantitative trait loci (eQTL) mapping between the genotype variants and transcriptome of three tissues (longissimus dorsi muscle, backfat, and liver) in 120 cattle. RESULTS: We identified 1,580 association signals for 21 beef agronomic traits using GWAS. We then illuminated 854,498 cis-eQTLs for 6,017 genes and 46,970 trans-eQTLs for 1,903 genes in three tissues and built a synergistic network by integrating transcriptomics with agronomic traits. These cis-eQTLs were preferentially close to the transcription start site and enriched in functional regulatory regions. We observed an average of 43.5% improvement in cis-eQTL discovery using multi-tissue eQTL mapping. Fine-mapping analysis revealed that 111, 192, and 194 variants were most likely to be causative to regulate gene expression in backfat, liver, and muscle, respectively. The transcriptome-wide association studies identified 722 genes significantly associated with 11 agronomic traits. Via the colocalization and Mendelian randomization analyses, we found that eQTLs of several genes were associated with the GWAS signals of agronomic traits in three tissues, which included genes, such as NADSYN1, NDUFS3, LTF and KIFC2 in liver, GRAMD1C, TMTC2 and ZNF613 in backfat, as well as TIGAR, NDUFS3 and L3HYPDH in muscle that could serve as the candidate genes for economic traits. CONCLUSIONS: The extensive atlas of GWAS, eQTL, fine-mapping, and transcriptome-wide association studies aid in the suggestion of potentially functional variants and genes in cattle agronomic traits and will be an invaluable source for genomics and breeding in beef cattle.

MAK: a machine learning framework improved genomic prediction via multi-target ensemble regressor chains and automatic selection of assistant traits.

Incorporating the genotypic and phenotypic of the correlated traits into the multi-trait model can significantly improve the prediction accuracy of the target trait in animal and plant breeding, as well as human genetics. However, in most cases, the phenotypic information of the correlated and target trait of the individual to be evaluated was null simultaneously, particularly for the newborn. Therefore, we propose a machine learning framework, MAK, to improve the prediction accuracy of the target trait by constructing the multi-target ensemble regression chains and selecting the assistant trait automatically, which predicted the genomic estimated breeding values of the target trait using genotypic information only. The prediction ability of MAK was significantly more robust than the genomic best linear unbiased prediction, BayesB, BayesRR and the multi trait Bayesian method in the four real animal and plant datasets, and the computational efficiency of MAK was roughly 100 times faster than BayesB and BayesRR.

A 1.1 Mb duplication CNV on chromosome 17 contributes to skeletal muscle development in Boer goats.

The Boer goat is one of the top meat breeds in modern animal husbandry and has attracted widespread attention for its unique growth performance. However, the genetic basis of muscle development in the Boer goat remains obscure. In this study, we identified specific structural variants in the Boer goat based on genome-wide selection signals and analyzed the basis of the molecular heredity of related candidate genes in muscle development. A total of 9 959 autosomal copy number variations (CNVs) were identified through selection signal analysis in 127 goat genomes. Specifically, we confirmed that the highest signal CNV (HSV) was a chromosomal arrangement containing an approximately 1.11 Mb (CHIR17: 60062304-61171840 bp) duplicated fragment inserted in reverse orientation and a 5 362 bp deleted region (CHIR17:60145940-60151302 bp) with overlapping genes (e.g., ARHGAP10, NR3C2, EDNRA, PRMT9, and TMEM184C). The homozygous duplicated HSV genotype (+/+) was found in 96% of Boer goats but was not detected in Eurasian goats and was only detected in 4% of indigenous African goats. The expression network of three candidate genes ( ARHGAP10, NR3C2, and EDNRA) regulating dose transcription was constructed by RNA sequencing. Results indicated that these genes were involved in the proliferation and differentiation of skeletal muscle satellite cells (SMSCs) and their overexpression significantly increased the expression of SAA3. The HSV of the Boer goat contributed to superior skeletal muscle growth via the dose effects of overlapping genes.

Genomic prediction of carcass traits using different haplotype block partitioning methods in beef cattle.

Genomic prediction (GP) based on haplotype alleles can capture quantitative trait loci (QTL) effects and increase predictive ability because the haplotypes are expected to be in linkage disequilibrium (LD) with QTL. In this study, we constructed haploblocks using LD-based and the fixed number of single nucleotide polymorphisms (fixed-SNP) methods with Illumina BovineHD chip in beef cattle. To evaluate the performance of different haplotype block partitioning methods, we constructed haploblocks based on LD thresholds (from r 2 > 0.2 to r 2 > 0.8) and the number of fixed-SNPs (5, 10, 20). The performance of predictive methods for three carcass traits including liveweight (LW), dressing percentage (DP), and longissimus dorsi muscle weight (LDMW) was evaluated using three approaches (GBLUP and BayesB model based on the SNP, GHBLUP, and BayesBH models based on the haploblock, and GHBLUP+GBLUP and BayesBH+BayesB models based on the combined haploblock and the nonblocked SNPs, which were located between blocks). In this study, we found the accuracies of LD-based and fixed-SNP haplotype Bayesian methods outperformed the Bayesian models (up to 8.54 ± 7.44% and 5.74 ± 2.95%, respectively). GHBLUP showed a high improvement (up to 11.29 ± 9.87%) compared with GBLUP. The Bayesian models have higher accuracies than BLUP models in most scenarios. The average computing time of the BayesBH+BayesB model can reduce by 29.3% compared with the BayesB model. The prediction accuracies using the LD-based haplotype method showed higher improvements than the fixed-SNP haplotype method. In addition, to avoid the influence of rare haplotypes generated from haplotype construction, we compared the performance of GP by filtering four types of minor haplotype allele frequency (MHAF) (0.01, 0.025, 0.05, and 0.1) under different conditions (LD levels were set at r 2 > 0.3, and the fixed number of SNPs was 5). We found the optimal MHAF threshold for LW was 0.01, and the optimal MHAF threshold for DP and LDMW was 0.025.

Cis-eQTL Analysis and Functional Validation of Candidate Genes for Carcass Yield Traits in Beef Cattle.

Carcass yield traits are of considerable economic importance for farm animals, which act as a major contributor to the world's food supply. Genome-wide association studies (GWASs) have identified many genetic variants associated with carcass yield traits in beef cattle. However, their functions are not effectively illustrated. In this study, we performed an integrative analysis of gene-based GWAS with expression quantitative trait locus (eQTL) analysis to detect candidate genes for carcass yield traits and validate their effects on bovine skeletal muscle satellite cells (BSCs). The gene-based GWAS and cis-eQTL analysis revealed 1780 GWAS and 1538 cis-expression genes. Among them, we identified 153 shared genes that may play important roles in carcass yield traits. Notably, the identified cis-eQTLs of PON3 and PRIM2 were significantly (p < 0.001) enriched in previous GWAS loci for carcass traits. Furthermore, overexpression of PON3 and PRIM2 promoted the BSCs' proliferation, increased the expression of MYOD and downregulated the expression of MYOG, which indicated that these genes may inhibit myogenic differentiation. In contrast, PON3 and PRIM2 were significantly downregulated during the differentiation of BSCs. These findings suggested that PON3 and PRIM2 may promote the proliferation of BSCs and inhibit them in the pre-differentiation stage. Our results further contribute to the understanding of the molecular mechanisms of carcass yield traits in beef cattle.

Improving Genomic Prediction with Machine Learning Incorporating TPE for Hyperparameters Optimization.

Depending on excellent prediction ability, machine learning has been considered the most powerful implement to analyze high-throughput sequencing genome data. However, the sophisticated process of tuning hyperparameters tremendously impedes the wider application of machine learning in animal and plant breeding programs. Therefore, we integrated an automatic tuning hyperparameters algorithm, tree-structured Parzen estimator (TPE), with machine learning to simplify the process of using machine learning for genomic prediction. In this study, we applied TPE to optimize the hyperparameters of Kernel ridge regression (KRR) and support vector regression (SVR). To evaluate the performance of TPE, we compared the prediction accuracy of KRR-TPE and SVR-TPE with the genomic best linear unbiased prediction (GBLUP) and KRR-RS, KRR-Grid, SVR-RS, and SVR-Grid, which tuned the hyperparameters of KRR and SVR by using random search (RS) and grid search (Gird) in a simulation dataset and the real datasets. The results indicated that KRR-TPE achieved the most powerful prediction ability considering all populations and was the most convenient. Especially for the Chinese Simmental beef cattle and Loblolly pine populations, the prediction accuracy of KRR-TPE had an 8.73% and 6.08% average improvement compared with GBLUP, respectively. Our study will greatly promote the application of machine learning in GP and further accelerate breeding progress.

Incorporating kernelized multi-omics data improves the accuracy of genomic prediction.

BACKGROUND: Genomic selection (GS) has revolutionized animal and plant breeding after the first implementation via early selection before measuring phenotypes. Besides genome, transcriptome and metabolome information are increasingly considered new sources for GS. Difficulties in building the model with multi-omics data for GS and the limit of specimen availability have both delayed the progress of investigating multi-omics. RESULTS: We utilized the Cosine kernel to map genomic and transcriptomic data as [Formula: see text] symmetric matrix (G matrix and T matrix), combined with the best linear unbiased prediction (BLUP) for GS. Here, we defined five kernel-based prediction models: genomic BLUP (GBLUP), transcriptome-BLUP (TBLUP), multi-omics BLUP (MBLUP, [Formula: see text]), multi-omics single-step BLUP (mssBLUP), and weighted multi-omics single-step BLUP (wmssBLUP) to integrate transcribed individuals and genotyped resource population. The predictive accuracy evaluations in four traits of the Chinese Simmental beef cattle population showed that (1) MBLUP was far preferred to GBLUP (ratio = 1.0), (2) the prediction accuracy of wmssBLUP and mssBLUP had 4.18% and 3.37% average improvement over GBLUP, (3) We also found the accuracy of wmssBLUP increased with the growing proportion of transcribed cattle in the whole resource population. CONCLUSIONS: We concluded that the inclusion of transcriptome data in GS had the potential to improve accuracy. Moreover, wmssBLUP is accepted to be a promising alternative for the present situation in which plenty of individuals are genotyped when fewer are transcribed.

Hereditary Basis of Coat Color and Excellent Feed Conversion Rate of Red Angus Cattle by Next-Generation Sequencing Data.

Angus cattle have made remarkable contributions to the livestock industry worldwide as a commercial meat-type breed. Some evidence supported that Angus cattle with different coat colors have different feed-to-meat ratios, and the genetic basis of their coat color is inconclusive. Here, genome-wide association study was performed to investigate the genetic divergence of black and red Angus cattle with 63 public genome sequencing data. General linear model analysis was used to identify genomic regions with potential candidate variant/genes that contribute to coat color and feed conversion rate. Results showed that six single nucleotide polymorphisms (SNPs) and two insertion−deletions, which were annotated in five genes (ZCCHC14, ANKRD11, FANCA, MC1R, and LOC532875 [AFG3-like protein 1]), considerably diverged between black and red Angus cattle. The strongest associated loci, namely, missense mutation CHIR18_14705671 (c.296T > C) and frameshift mutation CHIR18_12999497 (c.310G>-), were located in MC1R. Three consecutive strongly associated SNPs were also identified and located in FANCA, which is widely involved in the Fanconi anemia pathway. Several SNPs of highly associated SNPs was notably enriched in ZCCHC14 and ANKRD11, which are related to myofiber growth and muscle development. This study provides a basis for the use of potential genetic markers to be used in future breeding programs to improve cattle selection in terms of coat color and meat phenotype. This study is also helpful to understand the hereditary basis of different coat colors and meat phenotypes. However, the putative candidate genes or markers identified in this study require further investigation to confirm their phenotypic causality and potential effective genetic relationships.

Integrating genomics and transcriptomics to identify candidate genes for subcutaneous fat deposition in beef cattle.

Fat deposition is a complex economic trait regulated by polygenic genetic basis and environmental factors. Therefore, integrating multi-omics data to uncover its internal regulatory mechanism has attracted extensive attention. Here, we performed genomics and transcriptomics analysis to detect candidates affecting subcutaneous fat (SCF) deposition in beef cattle. The association of 770K SNPs with the backfat thickness captured nine significant SNPs within or near 11 genes. Additionally, 13 overlapping genes regarding fat deposition were determined via the analysis of differentially expressed genes and weighted gene co-expression network analysis (WGCNA). We then calculated the correlations of these genes with BFT and constructed their interaction network. Finally, seven biomarkers including ACACA, SCD, FASN, ACOX1, ELOVL5, HACD2, and HSD17B12 were screened. Notably, ACACA, identified by the integration of genomics and transcriptomics, was more likely to exert profound effects on SCF deposition. These findings provided novel insights into the regulation mechanism underlying bovine fat accumulation.

Identification of Candidate Genes Regulating Carcass Depth and Hind Leg Circumference in Simmental Beef Cattle Using Illumina Bovine Beadchip and Next-Generation Sequencing Analyses.

Genome-wide association studies are a robust means of identifying candidate genes that regulate economically important traits in farm animals. The aim of this study is to identify single-nucleotide polymorphisms (SNPs) and candidate genes potentially related to carcass depth and hind leg circumference in Simmental beef cattle. We performed Illumina Bovine HD Beadchip (~670 k SNPs) and next-generation sequencing (~12 million imputed SNPs) analyses of data from 1252 beef cattle, to which we applied a linear mixed model. Using a statistical threshold (p = 0.05/number of SNPs identified) and adopting a false discovery rate (FDR), we identified many putative SNPs on different bovine chromosomes. We identified 12 candidate genes potentially annotated with the markers identified, including CDKAL1 and E2F3, related to myogenesis and skeletal muscle development. The identification of such genes in Simmental beef cattle will help breeders to understand and improve related traits, such as meat yield.

Comparative transcriptomic analysis reveals region-specific expression patterns in different beef cuts.

BACKGROUND: Beef cuts in different regions of the carcass have different meat quality due to their distinct physiological function. The objective of this study was to characterize the region-specific expression differences using comparative transcriptomics analysis among five representative beef cuts (tenderloin, longissimus lumborum, rump, neck, chuck). RESULTS: We obtained 15,701 expressed genes in 30 muscle samples across five regions from carcass meat. We identified a total of 80 region-specific genes (RSGs), ranging from three (identified in the rump cut) to thirty (identified in the longissimus lumborum cut), and detected 25 transcription factors (TFs) for RSGs. Using a co-expression network analysis, we detected seven region-specific modules, including three positively correlated modules and four negatively correlated modules. We finally obtained 91 candidate genes related to meat quality, and the functional enrichment analyses showed that these genes were mainly involved in muscle fiber structure (e.g., TNNI1, TNNT1), fatty acids (e.g., SCD, LPL), amino acids (ALDH2, IVD, ACADS), ion channel binding (PHPT1, SNTA1, SUMO1, CNBP), protein processing (e.g., CDC37, GAPDH, NRBP1), as well as energy production and conversion (e.g., ATP8, COX8B, NDUFB6). Moreover, four candidate genes (ALDH2, CANX, IVD, PHPT1) were validated using RT-qPCR analyses which further supported our RNA-seq results. CONCLUSIONS: Our results provide valuable insights into understanding the transcriptome regulation of meat quality in different beef cuts, and these findings may further help to improve the selection for health-beneficial meat in beef cattle.

Transcriptional states and chromatin accessibility during bovine myoblasts proliferation and myogenic differentiation.

OBJECTIVES: Although major advances have been made in bovine epigenome study, the epigenetic basis for fetal skeletal muscle development still remains poorly understood. The aim is to recapitulated the time course of fetal skeletal muscle development in vitro, and explore the dynamic changes of chromatin accessibility and gene expression during bovine myoblasts proliferation and differentiation. METHODS: PDGFR- cells were isolated from bovine fetal skeletal muscle, then cultured and induced myogenic differentiation in vitro in a time-course study (P, D0, D2,and D4). The assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were performed. RESULTS: Among the enriched transcriptional factors with high variability, we determined the effects of MAFF, ZNF384, and KLF6 in myogenesis using RNA interference (RNAi). In addition, we identified both stage-specific genes and chromatin accessibility regions to reveal the sequential order of gene expression, transcriptional regulatory, and signal pathways involved in bovine skeletal muscle development. Further investigation integrating chromatin accessibility and transcriptome data was conducted to explore cis-regulatory regions in line with gene expression. Moreover, we combined bovine GWAS results of growth traits with regulatory regions defined by chromatin accessibility, providing a suggestive means for a more precise annotation of genetic variants of bovine growth traits. CONCLUSION: Overall, these findings provide valuable information for understanding the stepwise regulatory mechanisms in skeletal muscle development and conducting beef cattle genetic improvement programs.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle.

BACKGROUND: A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle. RESULTS: Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits. CONCLUSIONS: We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Transcriptome Analysis of Bovine Rumen Tissue in Three Developmental Stages.

Rumen development is a crucial physiological challenge for ruminants. However, the molecular mechanism regulating rumen development has not been clearly elucidated. In this study, we investigated genes involved in rumen development in 13 rumen tissues from three developmental stages (birth, youth, and adult) using RNA sequencing. We identified that 6,048 genes were differentially expressed among three developmental stages. Using weighted correlation network analysis, we found that 12 modules were significantly associated with developmental stages. Functional annotation and protein-protein interaction (PPI) network analysis revealed that CCNB1, CCNB2, IGF1, IGF2, HMGCL, BDH1, ACAT1, HMGCS2, and CREBBP involved in rumen development. Integrated transcriptome with GWAS information of carcass weight (CW), stomach weight (SW), marbling score (MS), backfat thickness (BFT), ribeye area (REA), and lean meat weight (LMW), we found that upregulated DEGs (fold change 0∼1) in birth-youth comparison were significantly enriched with GWAS signals of MS, downregulated DEGs (fold change >3) were significantly enriched with GWAS signals of SW, and fold change 0∼1 up/downregulated DEGs in birth-adult comparison were significantly enriched with GWAS signals of CW, LMW, REA, and BFT. Furthermore, we found that GWAS signals for CW, LMW, and REA were enriched in turquoise module, and GWAS signals for CW was enriched in lightgreen module. Our study provides novel insights into the molecular mechanism underlying rumen development in cattle and highlights an integrative analysis for illustrating the genetic architecture of beef complex traits.